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Abstract 

We consider deformations of the Whitham systems in the case when the initial 
system is close to linear one. It appears that the almost linear case requires a special 
procedure of the deformation of the Whitham system to make all the constructions 
stable in the linear limit. We suggest here a special deformation scheme which gives 
a stable deformation of the Whitham system for an almost linear initial system. 

1 General deformation schemes. 

We will consider deformations of the Whitham systems for non-linear differential equations 
which are close to linear ones for small values of some "non-linearity parameters" A. 

As is well known the Whitham method ([U |2j [3]) is connected with the slow modula- 
tions of periodic or quasiperiodic m-phase solutions of nonlinear systems 

F*{<p,<p t ,<p x ,...) = , i = l,...,n, cp = (<f 1 ,...,<f n ) (1.1) 
which are represented usually in the form 

<f\x,t) = ¥(k(V)x + w(U)t + O , U) (1.2) 

In these notations the functions k(U) and w(U) play the role of the "wave numbers" 
and "frequencies" of m-phase solutions and 6q are the initial phase shifts. The parameters 
of the solutions U = (U 1 , . . . , U ) can be chosen in arbitrary way, however, we assume 
that they do not change under arbitrary shifts of the initial phases 6q of solutions. 

The functions <& l (0) satisfy the system 

F l (<&,u a <&e«,kP<f> e e,...) = , i = l,...,n (1.3) 

and we choose for every U some function <&(0, U) as having "zero initial phase shifts". A 
full set of m-phase solutions of (11 .ip can then be represented in the form (11.21) . For m-phase 
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solutions of $n$ we have then k(U) = (^(U), . . . , k m (U)), w(U) = (^(U), . . . , w m (U)), 
#o = • • • , # m ), where U = (U 1 , . . . , U N ) are parameters of the solution. We require 
also that all the functions U) are 27r-periodic with respect to every 9 a , a = 1, . . . , m. 

Let us denote by A the family of the functions 3>(#, U) which depend on the parameters 
U in a smooth way and satisfy system (II. 3p for all U. We will assume also that A is the 
maximal family having these properties. 

As it is well known the famous area of the theory of integrable systems connected with 
m-phase solutions was started by the fundamental work of S.P. Novikov ([H]) where the 
algebro-geometric approach to the theory of quasi-periodic solutions was invented. An 
introduction of the finite-gap potentials as the stationary points of higher KdV-equations 
and the investigation of their algebro-geometric properties has become a corner stone of 
the algebro-geometric approach in the theory of solitons. As is well known by now the 
theory of algebro-geometric solutions has become much larger area then the pure theory 
of integrable systems and there are many other wide areas of mathematical physics where 
the varieties of Novikov potentials play a basic role. The algebro-geometric methods have 
become very efficient also in the theory of slow modulations of quasi-periodic solutions 
and play in fact the fundamental role in their consideration. 

We are going to consider here dispersion corrections to the Whitham systems and pay 
the special attention to the "almost linear" case. 

In Whitham approach the parameters U become slow functions of x and t: U = 
U(X, T), where X = ex, T = et (e -> 0). 

The functions U(X, T) should satisfy in this case some system of differential equations 
(Whitham system) which makes possible the construction of the corresponding asymptotic 
solution. More precisely (see [I]), we try to find the asymptotic solutions 

^(0,X,T) = ^& w (?(^ + 0,X,t\ e k (1.4) 

fc>o ^ e ' 

(where all ^rje) are 27r-periodic in 0) which satisfy the system (11.11) . i.e. 

F l (cp,e<p T ,etpx,- ■ •) = , i = l,...,n 

The function S(X,T) = (S 1 (X, T),..., S m (X, T)) is called a "modulated phase" of 
solution fTQl) . 

It is easy to see that the function ^^(0, X, T) should belong to the family of m-phase 
solutions of ( 11. ip at every X and T. We have then 

tf (O )(0,X,T) = &(0 + (X,T),V(X,T)) (1.5) 

and 

S%(X,T) = u a {XJ) , S%(X,T) = k a {U) 

as follows from the substitution of (II. 4p into system (II. ip . 

The functions ^f^){0, X, T) are defined from the linear systems 
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L^ 0o] (X,T)^ k) (e,X,T) = r w {0,X,T) (1.6) 

where L 1 -,^ e ^{X,T) is a linear operator given by the linearization of system ( II. 31) on 
solution (11. 5p . The resolvability conditions of systems (II. 6p can be written as the or- 
thogonality conditions of the functions i^){0,X, T) to all the "left eigen vectors" (the 

eigen vectors of adjoint operator) %l xr )](S + Oq(X,T)) of the operator U.^j & ^{X,T) 
corresponding to zero eigen-values. The resolvability conditions of (jl.6p for k = 1 

L) [VM {X,T)V {1) {0,X,T) = r (1) (0,X,T) (1.7) 

together with the relations k% = give the Whitham system for m-phase solutions of 
(II. ip which plays the central role in the slow modulations approach. 

Let us say that the resolvability conditions of (jl.6p can in fact be rather complicated 
in a general multi-phase case. Indeed, we need to investigate the eigen-spaces of the 
operators £[u,0 o ] an d ^nj e ] on ^ ne s P ace °f 27r-periodic functions which can be rather 
non-trivial in the multi-phase situation. Thus even the dimensions of kernels of Lnj,e ] 
and L, v g , can depend in non-smooth way on the values of U so we can have a rather 
complicated picture on the U-space ([151 HH1 E21 133]). 

These difficulties do not usually appear in the one-phase situation (m = 1) where the 
behavior of eigen-values of £[u,0 o ] anc ^ -^[U0 O ] * s usuan y ra ther regular. It is natural to 
introduce some regularity conditions on the space of one-phase solutions in this situation 
which will play an important role in the construction of asymptotic series (II .4p . We will 
assume here that the parameters k and u can be considered (locally) as the independent 
parameters on the family A and the total family of solutions of (II. 3p depends on iV = 2 + s 
(s > 0) parameters U u and the initial phase 6*0. 

Easy to see then that the functions &e(0 + 9 (X,T),XJ(X,T)) and 

$>{6 + 9o(X,T),XJ(X,T)) where £ is any vector in the space of parameters U u tan- 
gential to the surface k = const, u = const belong to the kernel of the operator L\ XT ^.. 

Let us represent the space of parameters U in the form U = (k,u,n) where k is the 
wave number, u is the frequency of one-phase solution, and n = (n 1 , . . . ,n s ) are some 
additional parameters (if they exist). 

Definition 1.1. 

We call the family A a full regular family of one-phase solutions of U.l\) if 

1) The functions &o{9, k, u, n), <fr n ;(#, k, u, n) are linearly independent and give the 
full basis in the kernel of the operator L^™ koJIl y 

2) The operator L\ 6o ku;n ^ has exactly s + 1 linearly independent "left eigen vectors" 
depending on the parameters U in a smooth way and corresponding to zero eigen-values. 
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The Whitham system is defined in this regular situation by the orthogonality con- 
ditions of the discrepancy f(y(9,X,T) to the "left eigen-vectors" k^ ui1 ^(9 + 9q), q = 
1 . . . . . ,s — 1 . 

As is well known, in many examples the Whitham system gives restrictions on the 
parameters U U (X,T) of zero approximation ^Sf^(9,X,T) and leaves free the parameter 
#o {X, T) . This fact was formulated in [UJ as a general Lemma for the case of the full 
regular family of (m-phase) solutions A. In fact, it is also true in a more general (m-phase) 
situation even without the requirements of " regularity" of the family A. Let us prove here 
the corresponding Lemma. 

Lemma 1.1. 

The orthogonality conditions of all the "left eigen vectors" °/-^}[u0 o ] corresponding 
to zero eigen values to the discrepancy f(i)(0,X,T) give restrictions on the functions 
U l/ (X,T) only and do not involve the functions 9q(X,T), . . ., 9™(X,T). 

Proof. 

Let us represent the main term of expansion (11.41) in the form 

4,Ue,x, T) = v(e+ s{x - T)+ f o{x ' T) ,v(x,t)} 

Easy to see then that the part of f(i) containing the functions 9q(X,T) has the form 

d'F i d'F i 

where the notations &F l /d'u a and d'F 1 /d'k a mean that we don't consider the dependence 
of the functions on k and u> in (II .3p and differentiate F l only with respect to the 
"explicit" k a and u a in (II .3p and then put if = $>(6, k, u>, n). 
However, the full derivatives 

F s ($,w a * fla ,...) , -^-F 4 ($,w a * flQ ,...) 



doj a K ' ° ' ' ' dk a 
(including the differentiation of <fr with respect to u> and k) are identically zero on the 
family A according to (II. 3p . So we can write that 



i.e. 



/(I) - £}[U,0 O ] X ' T ) ^OT + ^}[U,fl ] ®i a X ' T ^ ®0X 

We have then that f(i)(0, X, T) always belongs to the image of ^}rue ] without any 
restriction on the functions 9q(X,T). 
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Lemma 1.1 is proved. 



It is easy to see also from the proof of the Lemma that the functions 0$(X, T) generate 
the additions 3^ 6q T and &k a @ox to the function ^f^(0,X, T) which in the main order 
is equivalent to the effective " renormalization" 

u a -)• u a + e6% T , k a -)• k a + e6% x 

of the parameters (w, k) of zero approximation $(# + 6 , U). 

According to the deformation procedure used in [17] the parameters (X,T) become 
in fact unnecessary after a "renormalization" of the phase S(X,T) when the unnecessary 
"renormalization freedom" disappears. 

The Whitham system is a so-called system of Hydrodynamic Type, which can be 
written in the form 

a;(U)u» = b;(v)u» (i.s) 

with some matrices -A(U) and -B(U). In generic case the system (11.81) can be resolved 
w.r.t. the time derivatives of U and written in the evolution form 

1% = V»(U)U£ , u = l,...,N (1.9) 

(where V = A^B). 

Lagrangian properties of the Whitham system were investigated by Whitham ([3]) 
who suggested also a method of " averaging" of a Lagrangian function to get a Lagrangian 
function for the Whitham system. 

Another important procedure is the procedure of "averaging" of local Hamiltonian 
structures suggested by B.A. Dubrovin and S.P. Novikov ([T8J [2H1 ED])- The Dubrovin 
- Novikov procedure gives a field-theoretical Hamiltonian structure of Hydrodynamic 
Type for system (11.91) with a Hamiltonian function having the hydrodynamic form H = 
J h(XJ)dX. The Dubrovin - Novikov bracket for system (I1.9P has the form 

{U U {X),U' M {Y)} = g v »(U)5'(X -Y) + b u ^(V)U^5(X -Y) (1.10) 

which is called also a local Poisson bracket of Hydrodynamic Type. 

The Hamiltonian properties of systems (II. 9p are strongly correlated with their inte- 
grability properties. Thus it was proved by S.P. Tsarev ([21]) that all the diagonalizable 
systems (11.91) having the Dubrovin - Novikov Hamiltonian structure can in fact be inte- 
grated (S.P. Novikov conjecture). Actually the same is true also for the diagonalizable 
systems (I1.9P having more general weakly-nonlocal Mokhov-Ferapontov or Ferapontov 
Hamiltonian structures. Let us say also here that the Dubrovin - Novikov procedure of 
averaging of local Poisson brackets can be generalized also to the weakly-nonlocal case. 

The construction of asymptotic series (ll.4p for the case of a full regular family of (one- 
phase) solutions of (II. ip can be represented in a regular way. Namely, provided that the 
Whitham system is satisfied we find the first correction ^f^(6,X, T) at every X and T 
from system (11.71) . The function ^f^(8,X, T) is defined modulo the linear combination 
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cW(X,T)* e {0,X,T) + J24 1] (X,T)^ nl (9,X,T) (1.11) 

1=1 

of the eigen- vectors of L l j, v e ^ corresponding to zero eigen- values. The coefficients c^ k '(X, T) 

and di(X, T) arising at every step k and the initial phase 6 Q (X, T) can be used to provide 
resolvability of systems ( 11. 6ft in the higher orders so we can hope to find recurrently all 
the corrections ^^)(d,X,T). 

The structure of the recurrent procedure can be constructed in regular way. Let us 
note first of all that it was pointed out by J.C. Luke that the values c^'(X, T) are actually 
not involved in the resolvability conditions of (II. 6p for k = 2 and the order k = 2 gives 
restrictions on the initial phase Qo(X, T) instead ([4], see also [15j EEl [271 EE])- 

Let us say that this statement can be generalized in fact for all the orders k > 1. Let 
us prove here the corresponding Lemma which will not actually require the "full regular 
family" A of solutions of ( II .ip and can be formulated in fact in the most general (m- 
phase) situation. We will assume just that all the corrections ^(i), ■ ■ -, *&(k) are found in 
asymptotic solution ( 11.41) in general m-phase case and we can write the general solution 
*(fc) of ( II. 6p in the form 

m 

* (fc3 (0, X, T) = *' (fc) (0, X,T) + J2 <& (X, T) * ga (0,X,T)+J2 4^ T) Q v (0, X, T) 

a=l V 

where ^'^(0, X, T) is normalized in some way. Here the functions $Q a (0, X, T) always 
belong to the kernel of the operator L[u,e ] and Qi'(0, X, T) denote all the other (linearly 
independent) vectors from Ker Lnj t g ]. In our notations we put here also 

*(0,X,T) = *{0 + {X,T),\J(X,T)) = *(o)(0,X,T) 
Lemma 1.2. 

The functions Ca^ (X, T) do not appear in resolvability conditions of system U. 6\) in 
the order k + 1 . 

Proof. 

Assume first that k > 2. Let us look at the terms in fk + i(d,X,T) which contain the 
functions c£\x, T). Let us divide these terms in three groups: f^ +1 , f^, fjS in the 
following way: 

I) There are terms corresponding to the "correction" of the value F l ((p, k a ^pe a 1 ^ i Pef ) i ■ ■ ■) 
as a result of corrections of Sf^o) in the fe-th order. To describe these terms it is convenient 
to use again the fact that the correction e k Ca\x, T) $>ffa(0, X, T) to *S?(p)(0,X,T) is 

equivalent (modulo the terms of order e 2k ) to the correction e h+1 Ca \X, T) to the phase 
S a (X,T) of * {o) (0,X,T). 
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So the corresponding corrections of k-th and k + 1-th orders to F l will have the form: 
a) 



/•2-7T p2tt xpi(a\ m A m QI 



in 



This correction has order e k but it does not appear in f( fc ) since it describes in fact the 
freedom of the choice of *S?/ k \(0, X, T) on the k-th step. It is equal to zero since all $>g a 
belong to the kernel of operator L[u,6» ]- 

b) In the same way as in Lemma 1.1 the corrections of order (k + 1) to F l can be 
represented in the form 

a=l a=l 

where the notations & mean again 

^ \<p=*(0,X,T) , ^ l v =*(e,x,r) 

We denote a correction of this type by — e fe+1 f^ +1 )(0, X, T) according to our notations. 
Again we can state that 

m m 

f i+l) = ^2 C aT ^ l j[V,0o] + ^2 C aX ^}[U,0 O ] 

a=l a=l 

in the same way as in Lemma 1.1 and we get then that f( 7 fe+1 ) always belongs to the image 
of operator Z/[u,6» ]- 

Let us consider now the other two groups of terms in f(k+i) containing Ca\x,T). 

II) The second group represent the correction of 
f(i)[\l/(o)](0, X, T) as a result of the corrections of SE^o) m the k-th order. The interesting 
part of this correction of order (A; + 1) has again the form 

Jo Jo £^J o) (0',X,T) a K } 1 ; (2tt)- 
and is equal in fact to e k+1 J2™=i c ( a\x,T) f{ 1)ea (0, X, T). So we have 

m 
a=l 

This correction has the form of a shift of the phase of fm and belongs to f(jfc+i)- 

III) The third group is generated by the terms in f(fc+i) which contain the functions 

T) and ^^{d,X,T). Easy to see that the interesting part of these terms can 
be written in the form 
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i.e. 



flint* xt\- r r 62Ft{o) i 

f {k+1) (0,X,T) --J o -J Q pij^l^(W) x 



d m #' d m 6" 



p2ir r2ir m 

f[ k %(0,X,T) = - / ... / $>W(X,T) < ta (0',X,T)x 



x 



WW) l *=* {e ' X ' T) ^ ' X ' T) (27T)- (27T)- 



where the distribution 0") gives an "integral representation" of the operator L 1 -^ &Q y 
In a translationally invariant (in 0) case it's not difficult to prove then the relation 



/•27T /-27T m Araan 



/0 



/o a=1 W 



m 



+ 

So we get that 

m 



a=l 



which belongs to the image of L 1 -. Thus we get the statement of the Lemma for k > 2. 

Let us consider now the case k — 1. Let us represent for simplicity the solution 
* (1) (0,X,T) in the form 

m 

* {1) (0,X,T) = *' {1) (0,X,T) + ^r, c W(X,T)*e°(0,X,T) 

a=l 

where the freedom connected with vectors Qi>(0,X,T) is included in &LJ0, X, T). 

It is not difficult to see then that in order e 2 we will have all the terms described above 
for k > 2 and one extra term: 

/ (2) (0,X,T) - - I ... I 5 ^(0')8^(0") l,fi= ^ e ' x ' T)X 



o Jo 



a,/3=l I J I J 



containing the functions c^'(X,T). 

Consider the expansion of the values F l (<S>, u a <&ea, . . .) under a small shift of phase 
6 —s- + 50 where SO = c* 1 ) 5z, 5z — > 0. First of all we know that all the orders of this 
expansion should be equal to zero on A since the values F J (<fr, u a $> qol , . . . ) rcixicLiii zero on 
A after this shift. 

The second order of this expansion consists of two parts. The first part is equal to 
— ffX 1 (5z) 2 while the second one is equal to 



i*2tv i*2~k 

'o Jo 



p2ir Xpi(Q\ m d m f)' 



We get then 

m 

/(T = L me ] E cU(X,T)cf\x,T)& eae ,(0,X,T) 

a, (3=1 

i.e. f£V belongs to ImL[xj,e ]- We obtain thus the statement of the Lemma for all k > 1. 

Lemma 1.2 is proved. 



For the case of a full regular family of (one-phase) solutions of system (II. ip we can 
see then a regular scheme of construction of asymptotic solution ( jl.4p if the Whitham 
system ( 11. 8ft is satisfied. Namely, we find the correction ^f^(9,X, T) from the system 
( I1.7P modulo linear combination (11. lip and then try to find the correction \f r ( 2 )(#, X, T). 
The resolvability conditions for system ( 11. 6p for k = 2 then give us restrictions on the 
functions 9 Q (X, T), d[ 1 \x, T), I = l,...,s. Provided that the corresponding conditions 
are satisfied we find the solutions ^^){9,X, T) modulo a linear combination of the same 
type. Now at every step k > 2 we will have the restrictions on the functions c^ fc_2 ^(X, T), 
d\ k (X, T) and obtain the solution ^f^){9, X, T) modulo the linear combination 

s 

c^(X,T)<S> e (9,X,T) + ^4 fc) (X,T)$ ni (^X,T) (1.12) 

i=i 

So we get a regular way of the recurrent construction of functions *&(k){9, X, T). 

The procedure described above looks quite natural. Indeed, the functions 9q(X, T), 
c^(X, T), c( 2 '(X, T), ... represent in fact corrections to the modulated phase S(X, T) 
while the functions d\ 1 \x, T), d\ 2 \x, T), ... represent corrections to the parameters 
n = (n 1 , . . . ,n s ) of the zero approximation ^f^(9,X,T). The Whitham system (11.81) or 
(I1.9P gives restrictions on the functions u(X,T) = St, k(X,T) = Sx and rii(X, T) as a 
resolvability condition of system (II. 7p for the first correction ^f^(9, X, T). It is natural 
then that the resolvability conditions of system (II. 6p for k > 2 give restrictions on the 
corrections 9 (X, T), c«(X,T), c^(X,T), ... and ^ (1) (X,T), d[ 2) (X } T) } ... to S(X,T) 
and n l (X,T) in a successive order. 
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The Whitham solution (11. 4p can be rewritten also in the form: 



X ,T,e) = ^( 5( ^ T ' £) + 9, S x (X,T,e), S T (X,T,e), n{X,T,e)J + 

E (^f^ + 0, X, t) e fc (1.13) 



+ 

fc>i 



where we allow the regular e- dependence 

S(X,T,e) = J]S {fc) (X,T)e fe , n l (X,T,e) = ^nJ fc) (X,T)e fc 

of the phase S and the parameters (k,u, n) of the zero approximation ^f^(9, X,T) such 
that 

k(X,T,e) = S x (X,T,e) , u(X,T,e) = S T (X,T,e) 

In this approach the functions ^^){9, X, T) can be normalized in some way while the 
functions S(k), n\ k ^ can be used to provide the resolvability conditions of systems (II .6p . 

For instance, the normalization of ^f^){9, X, T) used in |UJ required that the main term 
of (I1.13P gives "the best approximation" to solution (jl.4p such that the rest of series f 1 1 . 1 3 j) 
is orthogonal to the functions 

^f S(X,T,e) + ^ Sx ( X ,T,e), S T (X,T,e), n(X,T,e) 



and 

* n i( S{X,T,e) +9, S x {X,T,e), S T {X,T,e), n(X,T,e] 



at every e. 

The functions ^^(9, X, T) satisfy linear systems analogous to (ll.6p . i.e. 

Vls^n^^T) = f~t k) (9,X,T) (1.14) 

The functions f«.) are slightly different from f^) since a "part of e-dependence" is 
included now in the main term of fll . 13j) . 

Let us formulate here the Lemma ([17]) about the systems on the functions S^)(X, T), 
nl k \(X,T), (k > 1) arising in this approach. 

Lemma 1.3. 

In the case of a full regular family of (one-phase) solutions of U.l\) the functions 
S(k){X, T), ri/ k )(X,T) satisfy the linearized Whitham system on the functions S(o)(X,T), 
ii( ) (X,T) with an additional right-hand part depending on the functions S^, . . ., S^-i), 

n (0); • • ■; n (fc-l)- 
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In fact, it is not difficult to show that the functions c ( - fc_1 ^(X, T), d\ (X, T), (k > I) 
satisfy rather similar systems in this case. 

Let us say that we can give a definition of a full regular family of m-phase solutions 
of fll.ip also for the case m > 1 if we require the properties of Definition 1.1 for generic 
k = (A; 1 , . . . , k m ) and uj = (u 1 , . . . , u m ). Let us give here the corresponding definition. 

Definition 1.1'. 

We call the family A a full regular family of m-phase solutions of ( ti.il) if 

1) The functions <3?(9«(0, k, uj, n) ; & n i(0, k, uj, n) are linearly independent and give for 
generic k and uj the full basis in the kernel of the operator L l -, eo kwn p 

2) The operator L^ 0o k u n j has for generic k and uj exactly m + s linearly independent 
"left eigen vectors" 

depending on the parameters U in a smooth way and corresponding to zero eigen-values. 

All the constructions described above can be used also for the full regular family of 
m-phase solutions if we require that the orthogonality conditions 

■■J o K^ XtT)] .(0 + o (X,T)) ft k) (0,X,T) — = (1.15) 

give the necessary and sufficient conditions of resolvability of systems (11. 6p . This is a 
serious requirement and it can be shown ((T5J HEl |32} [33] ) that it is not satisfied in the 
general case. However, there exist examples where this requirement is satisfied in m-phase 
situation so the same procedure of construction of asymptotic series ( 11.411 (or (11. 13ft can 
be used in this case. 

Let us assume in our further considerations that we have a full regular family of one- 
phase or m-phase solutions of (II .ip and the compatibility conditions of system (II .6p are 
defined by orthogonality conditions (|1.15p . Let us say again, that this situation is given 
by rather specific examples in a multi-phase case. 

As far as we know dispersive corrections to the Whitham systems were first considered 
by M.Y. Ablowitz and D.J. Benney ([5], also [6]- [7]) where the first consideration of 
a multi-phase Whitham method was also made. As was pointed out in [5] the higher 
corrections in Whitham method satisfy more complicated equations including " dispersive 
terms" and the Whitham system (ll.8p should in fact contain also the higher derivatives 
("dispersion") being considered in the next orders of e. 

In [4"Tt HH] a general procedure of deformation of the Whitham systems based on a 
renormalization of parameters was suggested. The deformations of the Whitham systems 
appeared in [4"Tt HH] have the so-called Dubrovin-Zhang form and were considered in 
connection with B.A. Dubrovin problem of deformations of Frobenius manifolds. Let us 
say here some words about this deformation scheme. 
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The higher corrections to Topological Quantum Field theories require the deformations 
( |38[ HOI H2]) of the Hydrodynamic Type hierarchies (11 .9p having the form 

W T = VZ(U)U» + ^^ ) (U,U x ,...,U fcX )e fc - 1 (1.16) 

k>2 

where all vY k \ are smooth functions polynomial in the derivatives Ux, • • •, U^x and having 
degree k according to the following gradation rule: 

1) All the functions /(U) have degree 0; 

2) The derivatives U% x have degree k; 

3) The degree of the product of two functions having certain degrees is equal to the 
sum of their degrees. 

Deformation (11. 16ft of system (11. 9ft implies also the deformation of the corresponding 
(bi-)Hamiltonian structures (ll.lOj) 

{U V (X),U»(Y)} = {U»(X),U»(Y)} + 

k 

+ H^ 1 E 5 (T )s (U,Ux,...,U (fe _ s) x) 6M{X-Y) (1.17) 

k>2 s=0 

where all B^ s are polynomial w.r.t. derivatives Uj, . . ., \J^- S )x and have degree (k — s). 

We call deformations of form (ll.16p - p.17p deformations of Dubrovin-Zhang type. Let 
us say that form (I1.16I) - (I1.17I) is not the only possible form of deformation of the Whitham 
system. For instance, a Lorentz-invariant scheme for nonlinear Klein-Gordon equation was 
considered in |49j . 

However, as we will see, the deformation procedure used in jlTJ HS1 HH] is not very 
good for "almost linear systems" where a non-linearity is rather small. This instability 
is connected with the general instability of the Whitham approximation in the higher 
orders for the case of the small amplitude of oscillations which was pointed out by A.C. 
Newell in [20]. In the next chapter we will suggest a deformation scheme of the Whitham 
system for such "almost linear" systems which should describe the slow modulations of 
periodic (or quasiperiodic) solutions in this situation. We will use here the deformations 
of Dubrovin-Zhang form f 1 1 . 1 6 j) although other types of gradation rules are also possible 
as well. 

In general a deformation of the Whitham system can be (in new notations) described 
in the following way: 

We look for a solution of ( II. ip having the form 

<p(0,X,T) = <&(S(X,T) + 0,S x ,S T ,n) + ^ * (fc) (S(X, T) + 0, X, T) (1.18) 

fc>i 

where the functions ^f(k) are now local functionals of (k, uj, n) and there derivatives having 
gradation degree k. We omit now the parameter e although we put first S = S(X,T,e), 
n = n(X, T, e). Now the higher derivatives of (k, uj, n) play the role of small parameters 



12 



in the expansion according to chosen gradation rule. The Dubrovin-Zhang gradation rule 
implies the following simple definitions: 

1) The functions k a (X,T) = S X (X,T), u a (X,T) = S$(X,T), and n l (X,T) have 
degree 0; 

2) Every differentiation with respect to X adds 1 to the degree of a function; 

3) The degree of the product of two functions having certain degrees is equal to the 
sum of their degrees. 

The functions ^f^){0, X, T) are defined from the linear systems 

L\ Sx , ST rtj%)(°i X , T ) = ftk)(0,X,T) (1.19) 

where f l ^(0, X, T) is the discrepancy having gradation k according to the rules introduced 
above. 

The functions V&a.) are uniquely normalized by the conditions 

/2-7T p2ir n rl m f) 
...jf o Sx.Sr.n) *\ k) (0,X,T) ^ = (1.20) 

J o - J o E $ n.(0. S x ,S r ,n) *\ k) (e,X,T) = (1.21) 

k > 1, (a = 1, . . . , m, Z = 1, . . . , s), and are local functionals of (k, u), n) and there 
derivatives having gradation degree k. 

The "renormalized" modulated phase S(X,T) and parameters n(X,T) satisfy now 
the deformed Whitham system 

Stt = (T? k) (k,u>,n,k x ,uj x ,n x ,...) (1.22) 

k>l 

n T = *7(*)( k ' w ' n ' k *' "x, ni, •••) (1-23) 

k>l 

where cr?U, rf, k s are general polynomials in derivatives kx, u>x, nx, kxx, ^xx, n-xx, ■ ■ ■ 
(with coefficients depending on (k, u>,n)) having degree fc. 

The functions Tjm coincide with the right-hand part of the Whitham system 
(11. ip and the functions <T(k), f1(k) are defined by orthogonality conditions arising on the 
k-th order of (11.191) . Let us remind again that we work here with a full regular family 
of (m-phase) solutions of (11.11) . Besides that we imply that the orthogonality conditions 
of f(fc) to all "regular" left eigen- vectors introduced in Definition 1.1' are equivalent 
to the resolvability conditions of the system (11.191) . The deformation of the Whitham 
system can be rewritten also in parameters (k, u>,n) in obvious way 

fvrp — (jj -yr 
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u t = Yl a h ( k > w > n ' k *> n x , • . . ) (1.24) 

k>l 

U T = V U W ' n ' kx ' nx ' • • • ) 

fc>l 

2 Deformation procedure for almost linear systems. 

Let us say however that the regular procedure of deformation formulated above does not 
behave well in the case of "almost linear" systems (11.11) . This means that the procedure 
of deformation and the corresponding asymptotic series (11.181) do not have a good limit 
for systems (11.11) when 

F*(A, (p, (ft, (f X} ■■■) = F*((p, tp t , (p X7 ...) + Fft\, ip, (p t , <p x , . . .) 

where Fq are linear in (cp, <p t , <p x , . . .) and i^(A, <p, ip t , <p x , ...)—> when A — > 0. 

A reason for such a behavior is that the operator Lj(\) can now be represented in the 
form 

L)(X) = % + %{\) 

where the operator -^o[ue ]i nas * n ^ ac ^ a ^ ar § er number of ("left" and "right") eigen- 
vectors corresponding to zero eigen- values than the operator I/| U6) u(A) (see [20], Chp. 
2). 

Indeed, let us consider for instance the system 

<Ptt ~ Vxx + ¥ + A</? 3 = , A -> (2.1) 
It is well known that system (12. ip has a two-parametric family of one-phase solutions 

<p = <&\{kx + Lot + 0q, /i) , (/i = u 2 — k 2 ) 

depending on the parameters k and u. 

Since the amplitude of the solution depends on /i in a singular way near the point A = 
it is more convenient to use the parameters k and A = & max — Q m in hi this situation. So 
we will write 

<p = $ x (kx + u(k,A,X)t +0 Q ,A) 

in the new notations. If we choose the initial phase such that the function <&\(9, A) has a 
local maximum at 9 = then we will have the additional condition $>\(9, A) = &\(—9, A). 
We assume also that <&\(9, A) = &\(9 + 2n, A) as usually and we have in the limit A — > 



lim $x(0,A) = A cos 9 
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The dependence n{A) disappears in the limit A — > and we have // = 1 for A = 0. 
The operator 

L(A) = ( u *-k?)^L + l + 3\& x (0) 

has only one (both "left" and "right") eigen-vector k\(9,A) = §\^(9,A) corresponding 
to zero eigen-value. We have 

lim $ xe (9,A) = -A sin 9 

which is an eigen-vector of the operator L . However, the function cos 9 is also an eigen- 
vector (both "left" and "right") of the operator L = d 2 /d9 2 + 1 corresponding to zero 
eigen-value. As a result, the operator L(\) has an eigen-vector (both "left" and "right") 
Cx(9, A) corresponding to a "small" eigen-value z/(A, A), such that is(\, A) — > for A — > 0. 

We obtain then that even if the resolvability conditions of systems (11.191) are satisfied 
the solutions ^/^)(9, X, T) can be singular at A — > if the right-hand parts f^)(9,X,T) 
are not orthogonal to the eigen-vector (\(9,A). We can see then that it's natural to 
require orthogonality of all f^) to both the eigen-vectors K\{9,A) and (\(9,A) to make 
the procedure regular in the limit A — > 0. 

In this situation we can not require anymore the normalization conditions f l 1 . 2 [) and 
just require that all the corrections ^/^) in expansion fl 1 . 1 8 [) are regular functions at 
A — > 0. Using these requirements we can obtain now both the deformed Whitham system 
for system (jl.ip and normalization conditions for coefficients c^ k \X,T) arising in the 
definition of the function ^^)(9, X, T). 

Let us consider now the corresponding procedure in a more general formulation. We 
will omit for simplicity the additional parameters n and consider a one-phase situation. 

Let us assume that we have an " almost linear" system (11. ip and a full regular family 
of (one-phase) periodic solutions 

ip\x,t) = $\(kx + u(k,A,\)t +0 Q ,A) 

Here we choose the parameters (k, A) instead of (k, u) in the limit A — > where A is 
some parameter playing the role of amplitude. We will also assume that the dependence 
oj(k,A, A) becomes the dispersion relation oj(°\k) for some branch of the spectrum of 
linear system 

Fi{<p,(p t ,<p x ,...) = (2.2) 
and the function <&\(9, k, A) becomes the corresponding solution A <&o(9, k) of the system 

F*(<f>,u(°\k)<£> e ,k<f> e ,...) = (2.3) 

It is convenient to assume also that all the systems (II. ip . (I2.2p . (12. 3 p are written in 
a real form and both the dispersion relation co^°\k) and the functions & l (9,k) are real 
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functions. We also have that 3>o(# + 2iT,k) = & (9,k) and the function $ ,ee(^, k) is 
proportional to $o($! k) in this case. 

The function &\ t o(9, k, A) is a "right" eigen-vector of the operator L(X) corresponding 
to zero eigen- value and we require that there is also exactly one "left" eigen vector K\\ t k,A] 
of L(X) corresponding to zero eigen-value for all A. The limit of $\ } e(9,k,A) for A 
is equal to A<&o,6»(#, k) and we denote Ko[k](9) the limit of the (normalized) vector 
K[\ t k,A}{@) for A — > 0. Besides that, we assume that there exist "left" and "right" real eigen- 
vectors C[A,fc,A](6 l ) and £[\,k,A] ($) °f the operator L(X) corresponding to a "small" eigen- 
value z/(A, k, A) which give the additional "left" and "right" real eigen- vectors Co[fc](#) an d 
£o[fc]($) °f operator L Q corresponding to zero eigen-value. As we said already the existence 
of the vectors €[\,k,A] an d C[\,k,A] is connected with the existence of additional real eigen- 
vectors $o,0fi/(0! k) (or <$>o(9, k)) and Ko[fc],6»(0) lying in the kernel of the operators L and 
Lq respectively. 

We assume that the branch (0) of the spectrum of L is non-degenerate and 

d'L (u,k) . n , _ 

— ^-^L=c(o) (fc) M ,k) ^ o 

(the notation d'/d'u means here that we don't keep u and k connected by the dispersion 
relation and consider them as free parameters after the substitution d/dt — > ud/d9, 
d/dx -> kd/dd). 

We have the identity 

Indeed, expression (|2.4j) is a limit of the expression 

2n (a s d'F i (X,<p,uipe,kipe,...) i d6 

^\[k,A{^,k)]i\ d ) Q^j \<p=4> x (0,k,A(u,k)) 2~ 

(in parameters (cj, k)) at A — > 0. This expression according to (11.31) is equal to 

27r d d6 

^\[k,A^,k)]i( e ) L)(\u } k) — & x (d,h,A{u,k)) — 

which is identically zero. 

However, for the vector £ [fc] we don't have an identity like (12.41) in generic situation, 
so we will imply 

Co[*]*(*) 3 £J - Lw ( » *5(M)f£ ^ (2.5) 
and 



o 



r-2n 



Cx[kM e ) gi^ W=^ x (e,k,A) 7^ = , A^O (2.6) 
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We can suggest now a recurrent procedure of construction of asymptotic solutions 
( I1.18P and the deformation of the Whitham system in the almost linear case. As we 
said already we are going to require the orthogonality of the discrepancies f(k)(0, X, T) to 
both the "left" eigen- vectors K X [k,A](9 + #o) an d Cx[k,A](@ + #o) & t eac h step. We don't put 
normalization conditions (ll.2(jp now and use the freedom in the coefficients c k (X,T) and 
9 (X, T) to provide orthogonality of fr k \ and Cx[k,A]- So we try to find a deformation of 
the Whitham system in the form 

k T = (u)(k,A)) x , A T = ^2a k (k,A,k x ,A x ,...) (2.7) 

k>l 

where all a k are polynomial in (k x , A x , k xx , A xx , . . .) and have degree k. Every function 
a k (k, A, k x , A x , . . .) is found as previously from the orthogonality conditions of ft^ and 
K\(9,X,T) in the k-th order. The system 

k T = (u(k,A)) x , At = a 1 (k,A,k x ,A x ) 

coincides with the Whitham system for the system (11.11) . 

The functions 6 (X,T) and c«(X,T), C ( 2 )(X,T), . . . are defined from the orthogonal- 
ity conditions of f(i) and f( 2 ), f(2), ... to the "left" eigen- vector £ A (6>,X, T) respectively. 
It's not difficult to obtain the form of the systems arising on the function 6q(X,T) and 
c^(X,T). Indeed, as we saw in the proof of Lemma 1.1 the part of f(i) containing the 
function 8o{X, T) has the form 

d'F i d'F 1 

f {1) (e,x,T) = - — 9 0T --^e 0X 

So we obtain that the equation on 6o{X, T) is the first order linear differential equation 
which can be written in the form 

Qx(k, A) 9 0T + P x {k, A) 9 0X = R x {k, A, k x , A x ) (2.8) 

where 

/" 27r f)' fPi Af) /" 27r fj' ipi Af) 

Q x (k,A) = J ( X[k>A]i (0)—(e,k,A)— , P x {k,A) = J ^^ — {B^A) — 

It can be also seen from the proof of the Lemma 1.2 that all the systems on the 
functions c (fc) (X,T) can be written in a similar form, i.e. 

Qx(k, A) c$ + P x (k, A) 45 + H x (k, A, k x , A x ) = R { x k) [k, A, c« c^} , k > 2 

(2.9) 

Q x (k,A)c$ + Px(KA)c^ x + H x (k,A,h x ,A x )cW + W x (k,A)c^c^ = R x 1] [k,A} 
where Q x (k, A) is a non- vanishing function at A — > according to (12. 6p . 
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For deformed Whitham system (12 .7p and a corresponding solution of the Cauchy 
problem it is natural to consider also the Cauchy problem for systems (I2.8l) - (l2.9p and 
define the functions 9q(X, T), c^ k '(X,T). So, finally we get a recurrent procedure of 
deformation of the Whitham system and the construction of asymptotic solutions (11.181) 
which is regular at A — > 0. Let us note that the deformed Whitham system here is 
different in general case from that described in the first chapter because of the different 
normalization of the corrections ^(fc) (6*, X, T). 

It is not difficult to understand why normalization conditions f l 1 . 2 [) are not very good 
in the almost linear situation. Indeed, in the case close to linear one the dependence 
A(k,u, A) becomes a singular function for A — > since k and u> are not independent 
in a linear case. So, the rigid fixation of phase S(X, T) by conditions (11.201) leads to 
quite an "unstable" behavior of the main term in (11.181) depending on St and Sx as on 
parameters. As a result, all the higher corrections in (11.181) become singular functions at 
A — > 0. Relations (I2.8p - (l2.9p describe in this case the corrections to S(X, T) which should 
not be included in the parameters of the main approximation to keep all the terms stable 
in the limit A — > 0. 

To avoid this difficulty arising in the Whitham approach it was suggested by A.C. 
Newell not to keep a rigid dispersion relation between u, k, and A and correct the 
Whitham equations in the almost linear case ( [20J ) . 

Indeed, we can see now that the most natural way is in fact not to put a rigid connection 
between the phase S(X, T) and the parameters (u, k) of the main approximation 

* {0) (9,X,T) = *(9 + 9 (X,T),k(X,T),A(k(X,T),cu(X,T))) 

if we don't want to introduce the additional parameters 9 (X, T), c^ k \X, T). 
Indeed we can try to write the asymptotic solution of (II. ip in the form 

p = $ {S{X, T) + e,k,A) + J2 *(*) (S(x, T ) + 0> X > T ) ( 2 - 10 ) 

k>l 

where we don't require the exact relations St = u>(k,A), Sx = k anymore. Instead, we 
allow now corrections to the phase S(X, T) and put these relations just in the main order 
of our asymptotic expansion. This approach gives us the possibility to put again the 
normalization conditions (11.201) and include the corrections generated by the functions 
9 (X,T), cW(X,T) to the phase S(X,T). 

The type of the corrections to the phase S(X, T) depends on a deformation type we 
choose. For a deformation of type (11.161) it is convenient to put 

S x = k , St = w(k,A) + ^2u {k) (k,A,k x ,A x ,...) (2.11) 

k>i 

where all uj(k) ar e polynomial in (kx, Ax, kxx, Axx, ■ ■ •) and have degree k. As previously, 
we assume here that any function of k and A has degree and every differentiation with 
respect to X adds 1 to the degree of a function. 
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Now we put again normalization conditions ( ll.20p and look for both the deformation 
of dispersion relation ( 12.1 ip and the Whitham system which has now the form 



k T = — lu(k,A) + ^2u ik) (k,A,k x ,A x ,- ■ ■)) , A T = ^a (fc) (k, A, k x , A x , ■ ■ •) 

\ k>l J k>l 

(2-12) 

System (I2.12p is a closed system on the functions k(X,T), A(X,T) and it's natural 
to consider this system as a full deformation of the Whitham system in the almost linear 
case. 

Now the procedure of construction of asymptotic series (I2.10p looks similar to the previ- 
ous one and we look for the functions ^ f (fc)(6 ) , X, T) as for local functionals of (k, A, k x , A x , . . .) 
polynomial in derivatives of (k,A) and having degree k. The functions *&(fy(9,X,T) sat- 
isfy systems analogous to ( I1.19P and are uniquely determined by normalization conditions 
f ll.20p . A difference in this approach is that we require now the orthogonality of all 
the functions f^)(9,X, T) to both the "left" eigen-vectors k,\(9,X, T) and £\(Q,X, T) 
of the operator Lj(X) corresponding to zero and "small" eigen- values. The functions 
u)(k)(k, A, k x , A x , . . .) and a^(k, A, k x , A x , . . .) are uniquely determined by the orthogo- 
nality conditions in the k-th order which gives the recurrent procedure of construction of 
the deformation of Whitham system ( I2.12p and asymptotic series ( I2.10p . 

Let us call system ( 12. 121) the deformation of the Whitham system in an almost linear 
case. 

In the limit A — > system (" 12 . 1 2[) gives a system describing slow modulations of solu- 
tions of a purely linear system (jl.ip . All the functions &o(9), Ko(9), Co($)> f(fc)($) coincide 
in this case with their first Fourier harmonics and in fact only two linearly independent 
functions of 9 ($ (6 ) ), $ o e (6 1 )) arise at every X and T after the substitution of "3?o(#) m 
initial system (II. ip . The orthogonality of all f^) to both the functions k, (9), Co(#) gives 
in this case the relations fr^ = and we have *&(k) = in this situation. 

All the sums in ( 12. lip and (I2.12p contain just the finite number of terms in the linear 
case and the exact solution of ( 11.11) is given by the relation 

<p(6,X,T) = A{X,T)$(S{X,T) + 9,k{X,T)) 

where the functions S(X,T), A(X,T), k(X,T) satisfy the systems (l2TTTl) - fl2TT2l) . 

For the A-expansion of systems (I2.1ip -( l2.12p the Fourier expansion of the functions 
<&x(9,k,A), K\[k,A](9), C\[k,A](9) can be used. Indeed, the higher harmonics of functions 
$a, K\, Ca decrease usually as some power of A so we need just a finite number of Fourier 
harmonics at a given order of A. This approach can be connected then with the well 
known method of including of a non-linearity to the slowly modulated solutions of linear 
systems (see for instance [5]). Thus, the method of derivation of modulation equations 
for the parameters of the first Fourier harmonic of the solution (|3J, Chp. 15-16) can be 
considered as a calculation of the first A-correction in the system (I2.11l) - (l2.12p in this case. 

Finally, let us say just some words about multi-phase situation. As we told above, 
a multi-phase situation can be much more complicated for the construction of asymp- 
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totic solutions f l 1 . 1 8 j) or f )2.10p . So all our remarks should be addressed to some special 
cases when the orthogonality of the discrepancies f(k)(0, X, T) to the "regular" left eigen- 
vectors K^L A j(0) of the operator L % -{\) is sufficient for the determination of the corrections 
*( k )(0,X,T). 

We will assume according to our definition that we have m linearly independent "reg- 
ular" left eigen-vectors k^L a i(0) of the operator L(X) corresponding to zero eigen-values 
which is equal to the number of functions &\ t ga(6, k, A) giving the "regular" right eigen- 
vectors of L(X) corresponding to zero eigen-values. We assume also that we have m 
parameters A = (A 1 , . . . , A m ) playing the role of amplitudes (say the amplitudes of the 
main Fourier harmonics {cos(6 a + 0$)} in the expansion of <&\(0)) and we can express 
the function u a in the form u a = oo a (k, A, A) for rather small A. It is natural to accept 
then that in the limit A — > we have 

m 
a=l 

where g a = (g", . . . , gffi are some constant vectors and the functions u/*(k, A, A) become 
m independent dispersion relations u a = u a (k a ) also independent of A. 

According to our general approach we will assume that we have m additional " regular" 
left real eigen-vectors d|kA](^)' (M — 1 ; ■ • ■ j^t.) of -^(A) corresponding to "small" eigen- 
values z/( g )(A, k, A) which give m additional left eigen-vectors Copki(^) °f the operator Lq 
corresponding to zero eigen-values. The number of vectors Q?A corresponds in this case 

to m additional right eigen-vectors g Q cos(# a + 9q) of the operator Lq corresponding to 
zero eigen-values. 

As in the one-phase situation we will require now orthogonality of the functions 
f(fc)(0,X, T) to both the sets of the left eigen-vectors {4]1a](^} anc ^ {CxrkAiW} anc ^ 
try to find a deformation of the Whitham system and the dispersion relation in the form 



La 

T ~ dX 



S a x = k a , S% = w«(k,A) + 5>^(k,A,k x ,A x ,...) (2.13) 

k>\ 

L a (k,A)+^^ ) (k,A,k^A x ,...)j , ^ = ^a^(k,A,k x ,A x ,...) 

V k>l J k>l 

(2.14) 

where all the functions a%. satisfy the same requirements as in the one-phase situa- 
tion. 

Using our assumptions we can try then to find the functions vI/oa as local functionals 
of (k, A) and their derivatives and repeat all the steps of one-phase situation. All the 
functions afy, uj? k s will be uniquely determined in this case. However, we should remind 
again that this assumption is rather serious in the m-phase situation and in fact is not 
valid in general case. 
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